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Abstract 

Bose-Einstein Correlations (BEC) of three identical charged pions were studied in 4 x 10^ 
hadronic Z" decays recorded with the OPAL detector at LEP. The genuine three-pion correlations, 
corrected for the Coulomb effect, were separated from the known two-pion correlations by a 
new subtraction procedure. A significant genuine three-pion BEC enhancement near threshold 
was observed having an emitter source radius of rj, = 0.580 ± 0.004 (stat.) it 0.029 (syst.) fm and 
a strength of = 0.504 it 0.010 (stat.) it 0.041 (syst.). The Coulomb correction was found to 
increase the value by ~ 9% and to reduce rj, by ~ 6%. The measured corresponds to a value 
of 0.707 ±0.014 (stat.) ±0.078 (syst.) when one takes into account the three-pion sample purity. 
A relation between the two-pion and the three-pion source parameters is discussed. 
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1 Introduction 

Bose-Einstein Correlations (BEC) between pairs of identical bosons, mainly the n^n^ system, have 
been extensively studied in a large variety of interactions and over a wide range of energies [[l|, ^. 
These correlations, which are present when the bosons are near to one another in phase space, are 
used to estimate the size of the emitter of the particles and more recently in QCD-based models to 
describe fragmentation and hadronisation in high energy reactions The two-pion BEC effect 

has lately also been discussed in connection with the measurement of the W mass in the reaction 
e+e- W+W" hadrons at LEP2. 

In systems of more than two identical bosons, BEC are also expected to be present. These higher- 
order correlations may affect the multi-hadron production and are also of interest in intermittency 
studies [p. Detection of the "genuine" multi -boson BEC is complicated by the fact that they have 
to be isolated from the lower-order boson correlations and therefore require large data samples. In 
addition, systems of several identical charged bosons, placed nearby in phase space, are subject to 
a relatively large repulsive Coulomb interaction which may suppress the BEC effect. As a conse- 
quence, only relatively few higher order BEC studies of three and more charged pions have been 
reported [0-[T^]. In those studies, mainly due to lack of statistics, it was not possible to isolate and 
verify the genuine multi-boson BEC from the lower-order ones. Attempts have been made to infer 
from the measured BEC of the multi-boson systems the individual contribution of each of the higher 
order correlations by model dependent formulae. A significant genuine BEC signal of three identical 
charged pions, n^n^n^, has been reported in a hadron-proton interaction experiment [ [T2| ] and more 
recently in a LEP experiment [|T3]] where, however, the Coulomb effect was neglected. 

Here we report on a BEC study of the tt^tt'^tc'^ system carried out with a large sample of approx- 
imately 4 X 10^ hadronic Z*^ decays, recorded by the OPAL detector at the e+e^ LEP collider during 
the years 1991 to 1995. In this analysis, which takes into account the Coulomb effect, we have iso- 
lated the genuine three-pion BEC and estimated the size of the emitter. In Section 2 we introduce the 
extension of the two-boson BEC to the system of three identical bosons. Section 3 is devoted to the 
procedure used for the Coulomb correction of the three-pion BEC and in Section 4 the experimental 
details are given. In Section 5 we describe our method for the extraction of the genuine n^n^n^ 
BEC and present the results obtained from our analysis. The relations between the two-pion and the 
genuine three-pion BEC parameters are explored in Section 6. Finally, the summary and conclusions 
are presented in Section 7. 

2 The three-boson correlation function 

In describing the three-boson BEC we follow the approach which was also adopted, for example, 
in [^. The BEC of pairs of identical bosons can be formally expressed in terms of the normalised 
function 



where o is the total boson production cross section, pi(pi) and do/ dpi are the single-boson density in 
momentum space and the inclusive cross section, respectively. Similarly p2ipi,P2) and d^o/dpidp2 
are respectively the density of the two-boson system and its inclusive cross section. The product of the 
independent one-particle densities pi(;?i)pi(p2) is referred to as the reference density distribution, or 




(1) 
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reference sample, to which the measured correlations are compared. The inclusive two-boson density 
PiiPiiPi) can be written as: 

P2(P1,P2) = Pl(j9i)pi(p2) +^2(Pl,P2) , (2) 

where ^2(^1,^2) represents the two-body correlations. In the simple case of two identical bosons the 
normalised density function 7?2, defined in Eq. |T], already describes the genuine two-body correlations 
and has been referred to in previous BEC studies of OPAL [[T4|,p3|] as the C2 correlation function. 
Thus one has ^ 

C2=R2 = 1+K2{PUP2) , (3) 

where K2 {pi,P2) = K2{pi,P2)/[Pi{pi)Pi{P2)] is the normalised two-body correlation term. Since 
Bose-Einstein correlation are present when the bosons are close to one another in phase space, one 
natural choice is to study them as a function of the variable Q2 defined by 

Ql = qh = -{pi-P2f=Ml- V , 

which approaches zero as the identical bosons move closer in phase space. Here pi is the four- 
momentum vector of the iih. particle, n is the boson mass and is the invariant mass squared of the 
two-boson system. 



In the parametrisation proposed by Goldhaber et al. [|T6|], C2 has the form 

C2(e2) = l+^2e-ei^^ (4) 

where r2 estimates the size of the two-boson emitter which is taken to be of Gaussian shape. The 
strength of the BEC effect, frequently referred to as the chaoticity parameter, is measured by X2 
which varies in the range < ?i2 < 1 • 

The inclusive density of three bosons, p3(pi,;?2,i?3), includes the three independent boson mo- 
mentum spectra, the two-particle correlations K2 and the genuine three-particle correlations ^3, namely: 

p3{PUP2,P3) = pl{pi)pl{p2)pi{p3)+Y,PiiPi)K2{Pj,Pk)+K3{pi,p2,P3) , (5) 

(3) 

where the summation is taken over all the three possible permutations. The normalised inclusive 
three-body density, is then given by 

„ P3(P1,P2,P3) 1 , D , 2^ ^ ^ 

^3 = ^ . . . ^ . ^ = i+Rl,2+K3{Pl,P2,P3) ■ (6) 
Pl{Pl)Pl{P2)Pl{p3) 

Here 

L{3)PliPi)K2{Pi,Pk) 

Ki 2 = } — ^ } — ^ 7 — C'J 

Pl{Pl)Pl{P2)Pl{p3) 

represents a mixed three-boson system in which only two of them are correlated, and 

~ r \ K3{pi,p2,P3) 

^3 (/'i,/?2,;?3) = 7 — T— 7 — T, (8) 

Pl{pi)pl{p2)pl{P3) 

represents the three-boson correlation. In analogy with C2, one can define a correlation function C3, 
which measures the genuine three-boson correlation, by subtracting from R3 the term which contains 
the two-boson correlations contribution. Thus 

C3 = R3-Ri,2 =1+^3 {PUP2.P3) , (9) 
5 



which depends only on the genuine three-boson correlations. For the study of the three-boson corre- 
lations we use the variable Qt, which is defined as 

(3) 

where the summation is taken over all the three different /, j boson pairs and is the invariant mass 
squared of the three-boson system. From the definition of this three-boson variable it is clear that as 
Q2, approaches zero so do the three qtj values which eventually reach the region where the two-boson 
BEC enhancement is observed. 

The genuine three-pion correlation function C^^iQ'i) can be parametrised by the expression [^] 

C3(e3) = 1+2x3^-^3'-?, (10) 

where ?i3, which can vary within the limits < ?i3 < 1, measures the strength of the three-boson BEC 
effect and estimates the size of the three-boson emitter. The factor two which multiplies X3 arises 
from the presence of two possible diagrams with exchange of all the identical pions within a triplet. 
To extract from the data values for the strength and the emitter size r3, we modified Eq. |T^to read 

C3(e3) = K(l + 2V"2k?)(l+e<23) , (11) 

where K is a normalisation factor and the linear term ( 1 + £^3 ) accounts for the long range correlations 
arising from charge and energy conservation and phase space constraints. Higher-order ^3 terms for 
the long range correlations were found not to be needed in the present analysis. 



3 The Coulomb correction 

The observed BEC of identical charged bosons is suppressed by the Coulomb repulsive force. To 
account for this effect a correction to the measured BEC distribution is required. If C2{Q2) is the 
two-boson correlation in the presence of the Coulomb effect then it is related to the true C2^^{Q2) 
through the function G2{Q2), so that 

C2{Q2) = G2{Q2)Cr%Q2) . (12) 

In the case that the reference sample is a Monte Carlo generated data without the Coulomb effect, 
^2(22) can be expressed by the Gamow factor [[TTp: 

G2(e2) = 27ni/(e2'^^-l), (13) 

where r\ = (Xsm^i^H-i/Qi- Here ei and e2 are the charges, in positron units, of the two bosons having 
a mass of /j, and ttem is the fine-structure constant. Recently alternative Coulomb corrections [[rS,|l9p 



for the two-boson system have been proposed, which are based on non-Gaussian parametrisations. 
These could not be extended to the three-pion system and therefore were not used in our analysis. 

For a given three charged bosons system, with boson pairs having Q2 values of qi^2, ^1,3 and ^2,3, 
the Coulomb correction G3 can be approximated, in terms of the G2 function by [|^]: 

^3(23) = (G2(^l,2)G2(^2,3)G2(^l,3)) , (14) 
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where the average is taken over all experimentally accessible values of qij which satisfy the condition 

The function GT,{Q-i) can be evaluated through Eq. |14| or by using a slightly more precise for- 
mulation proposed in In our analysis the differences between the results of these two methods 
were smaller than the statistical errors so that the simpler method was sufficiently accurate for our 
purposes. 



4 Experimental setup and data selection 
4.1 The OPAL detector 

Details of the OPAL detector and its performance at the LEP e+e collider are given elsewhere [^. 
Here we will describe briefly only those detector components pertinent to the present analysis, namely 
the central tracking chambers. 

Besides a silicon microvertex detector, the central tracking chambers consist of a precision ver- 
tex detector, a large jet chamber, and additional z-chambers surrounding the jet chamber. The vertex 
detector is a 1 m long, two-layer cylindrical drift chamber that surrounds the beam pipeQ. The jet 
chamber has a length of 4 m and a diameter of 3.7 m. It is divided into 24 sectors in (|), each equipped 
with 159 sense wires parallel to the beam ensuring a large number of measured points even for parti- 
cles emerging from a secondary vertex. The jet chamber also provides a measurement of the specific 
energy loss, dE /dx, of charged particles A resolution of 3 — 4% on dE jdx has been obtained, 
allowing particle identification over a large momentum range. The z-chambers, 4 m long, 50 cm wide 
and 59 mm thick, allow a precise measurement of the z-coordinate of the charged tracks. They cover 
polar angles in the region | cos(0) | < 0.72 and 94% of the azimuthal angular range. All the chambers 
are contained in a solenoid providing an axial magnetic field of 0.435 T. The combination of these 
chambers leads to a momentum resolution of a^,/ pt ~ a/ (0.02)^ + (0.0015 • ptY, where pt in GeV/c 
is the transverse momentum with respect to the beam direction. The first term under the square root 
sign represents the contribution from multiple Coulomb scattering W%- 



4.2 Data selection 

The analysis was performed with the OPAL data collected at LEP with centre of mass energies on 
and around the peak with the requirement that the jet and z-chambers were fully operational. 
The hadronic decays were selected according to the number of charged tracks and the visible 
energy of the event [|24|]. We applied the same track quality and dE jdx cuts described in a former 
OPAL BEC study of two identical charged pions [[T4l]. Furthermore, events with a thrust angle of 
I cos(0thmst) I > 0.82 with respect to the beam axis were rejected. Finally we accepted only events 
with a relative charge balance of — < 0.25, where n+ and nT are respectively 

the observed numbers of positively and negatively charged tracks. Following these criteria a total of 
2.65 X 10^ hadronic decay events were used in the analysis. 

To avoid configurations with overlapping tracks, we rejected pion pairs if their invariant mass was 

'a right-handed coordinate system is adopted by OPAL, where the x-axis points to the centre of the LEP ring, and 
positive z is along the electron beam direction. The angles 9 and (|) are the polar and azimuthal angles respectively. 
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less than 0.4 GeV and if the opening angle between them in the plane perpendicular to the beam axis 
was less than 0.05 rad. In the present analysis all charged tracks are assumed to be pions. From Monte 
Carlo (MC) studies we have estimated that the average pion purity of our charged track sample is 
89.3% with a systematic uncertainty of ±2.2% and a negligible statistical error. Thus the pion purity 
of the three charged track system is71.3±5.3%. 



5 Analysis and results 

The BEC analysis used the hadronic TP decay data of OPAL where in each event all possible tt^tt^tc^ 
combinations were taken as the data sample. For the reference distribution we used a Monte Carlo 
generated sample p5| ] of 4 x 10^ JETSET 7.4 events which have passed a full detector simulation 
but do not include BEC and Coulomb effects. This JETSET Monte Carlo program, which includes 
most of the known resonances which decay into the tt+tt^ and tt^tt^tt^ final states, was carefully 
tuned to the OPAL data [ETp. Here one should note that the TT^TT'^^Tr^ systems of the measured data 



contain one pair of identical pions and therefore cannot be used as a reference sample. Thus: 



where /^^^^( is the number of the tt^tt^tt^ data combinations and M^^^ (23) is the correspond- 
ing number of Monte Carlo tt^tt^tt^ combinations at the same Qj, value. The Monte Carlo sample 
was normalised to the data in a region far away from any observable BEC enhancement. This 
was achieved by requiring that the integrated number of the Monte Carlo entries in the Qj, range 
1.6 — 2.0 GeV was equal to that of the data. The measured ^3(23) distribution is shown in Fig. [l]a 
in the range of 0.2 < Q3, < 2.0 GeV. The data points below 0.25 GeV have relatively large errors, 
and a lower three-track separation efficiency of identical charged tracks. Therefore a lower limit of 
0.25 GeV was imposed on the analysis. In the figure a clear enhancement is observed in the region 
below 23 = 1 GeV. This enhancement can be interpreted as coming from both the known two-pion 
and from the genuine three-pion BEC. 



5.1 The extraction of the genuine three-pion BEC 

To extract the genuine three-pion BEC one has to subtract fromi?3(23) the contribution coming from 
the well known two-pion BEC. In our analysis this last contribution is evaluated from the mixed- 
charged Tt^Tt^if^ combinations of the data. To this end, it is convenient to rewrite R-i{Q-i) as follows: 

R(n\ ^^^^(Q3) . , ^^^^(23)-M^^^(23) . , 5^^^(23) 

The total excess 5^^^ (23) above the Monte Carlo expectation has two contributions. The first from 
the two-pion BEC, §2 ^^(23), and the second from the genuine three-pion BEC, ^g^^i^^iOs), so that: 

5±±±(23) = §±^^(23) +5j±4,(23) . (17) 
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Figure 1: The n^n^TZ^ BEC measured distributions, without Coulomb correction, as a function of 
Q3. (a) the measured ^^3(23) and i?i,2(G3) distributions before subtraction of the two-pion BEC 
and (b) the C3(23) distribution after the subtraction of the two-pion BEC. These measurements are 
represented by points with statistical error bars. The solid line in (b) represents the fit result of Eq. [11], 
in the range 0.25 < Q3 < 2.0 GeV, to the measured C3(23) distribution. 
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From this it follows that what should be subtracted from RsiQs) to obtain the genuine three-pion 
BEC, is: 

Because we utilise for the subtraction the same events used for the ^3(23) measurement, there exists 
in every event with given charged multiplicity m and charge balance A = |n+ — | values a relation 
between 5j^^(23) and d^^^iQs) = N^^^Qs) -M^^^iQs). Here N^^^iQs) is the number of 
n^n^TT^ combinations in the data and M^^^{Qt,) is the corresponding number of combinations for 
the MC generated sample properly normalised to the data in the range of 1 .6 — 2.0 GeV. 

If we define n^^^ as the number of n^n^nr^ combinations in a given event, then it can be related 
to n^^"^, the number of n'^n^n^ combinations. A straightforward combinatorial calculation yields 
that for an event with given m and A values, one has: 

^++=3 1+4^- — ^ . (19) 

In our analysis we utilise, after the application of the appropriate cuts, all the hadronic decays 
lying within a wide range of multiplicity and charge balance. If we define {n^^^ /n^^^) as the ratio 
n^^^ /n^^^ averaged over all multiplicity and charge balance values in our data, then: 

j 6^^^{Q3)dQ3 = X / §±±=^(23)^23 . (20) 

The integrations are carried out over the Q3 region where 62 ^^(23) and 6^^^ (23) are different from 
zero. We further verified from MC studies that, to a good approximation, relation ( p0| ) holds also in 
its differential form, that is: 

5j^^(23) = (|[^^^x5±±T(23). (21) 

The average {n^^^ / n^^^) depends on the multiplicity m and the charge balance A distributions of 
the data sample. Using the MC hadronic decay sample we determined, by counting the number of 
combinations in the 23 range of 1.6 — 2.0 GeV, that {n^^^ / n^^"^) =3.69 with a negligible statistical 
error. The same value is obtained when the MC sample is replaced by the data sample. This value 
shifts to 3.71 when the 23 range is enlarged to 1.5 — 2.0 GeV. We also studied the variation of this 
ratio on the pion purity. To this end we evaluated this ratio from the MC sample using only tracks 
which were generated as pions with the result that (n^^^/n^^^) = 3.70. The effect of the shift from 
3.69 to 3.71 on the BEC parameters was found to be negligible in comparison to the statistical errors 
and to other systematic uncertainties (see Table ^. 

In Fig. 1^ we show the 5^^^(23) and the (n^^^/n^^^) x 5^^^(23) distributions as a function 
of 23- As can be seen, the two distributions are similar in the higher 23 region, namely between 0.7 
and 2.0 GeV. The slight difference between the two distributions can be attributed to the systematic 
uncertainties given in Table ^ in particular those listed as items (f) and (g). Thus in this 23 range 
the excess of ti^ti^tt^ combinations is fully accounted for by the excess seen in the n^n^n^ due to 
the two-pion BEC. In the lower 23 region an excess of 6^^^ (23) over (n^^^/n^^^) x 6^"^^ (23) 
is observed which can no longer be attributed to the two-pion BEC and is therefore identified as the 
genuine three-pion BEC contribution, Sgg^^j^g(23)- Thus: 



C3(23) = i?3(23)-i?l,2(23) = M±±±(23)~M±±±(23) \^/ ' ^ ^ 
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Figure 2: (a) 5^^^(23) and (b) (n^^^/n^^^) x 5^^''^(23) as functions of Q3. The errors plotted are 
the statistical ones. The difference between the distributions (a) and (b), in the range Qt, < 0.7 GeV, 
is due to the genuine three-pion BEC. 



In the present analysis we have used Eq. to subtract the contributions due to the two-pion BEC. 

The measured distribution i?i,2(23) is shown in Fig. [I]a and the resulting C^{Q3) distribution is 
shown in Fig. |T]b, where a significant genuine three-pion BEC enhancement is clearly present. The 
solid line in Fig. |I]b represents the fit result of Eq. [iT] to the data. The fitted values of X3 and rj, and 
the correlation factor px, ,- are given in Table |l] together with the %^ value divided by the number of 
degrees of freedom (d.o.f.). 



5.2 Evaluation of the Coulomb effect 

The Coulomb correction to the genuine BEC, defined as G3(23) in Eq. [T^, can be applied either to 
the data or to the MC reference sample. These two possibilities are not expected to yield identical 
results since, unlike the MC generated sample, the data are affected by both the BEC and the Coulomb 
interactions. In our BEC analysis we chose to apply the Coulomb correction to the data and utilised 
the results coming from the second possibility as a measure of the systematic errors. 

The Coulomb correction was accounted for by assigning to every three-pion combination of the 
data a weight equal to 1/ (G2(^i.2)G2(^2.3)G2(^i.3)). This automatically assures that 03(23), defined 
in Eq. [1^, is averaged only over all accessible experimental values of qtj. Using this procedure, the 
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Parameter 


Without Coulomb Corr. 


With Coulomb Corr. 


h 


0.462 ±0.012 


0.504 ±0.010 


rs [fm] 


0.616 ±0.005 


0.580 ±0.004 


K 


1.003 ±0.001 


1.026 ±0.001 


£ [GeV 1] 


-0.001 ±0.001 


-0.015 ±0.001 


PX.r 


+0.887 


+0.883 


xVd.o.f. 


218/171 


190/171 



Table 1: Results of the fit of Eq. [IT] to the measured C3(23) distributions without (Fig. |T]b) and with 
(Fig. ^) Coulomb correction carried out over the range of 0.25 < Q3 < 2.0 GeV. The fitted values 
are the genuine three-pion emitter and the BEC strength together with the long range correlation 
parameter £ and the normalisation factor K. The errors are the statistical ones obtained by the fits, and 
Px,r is the correlation factor between X3 and r^. The quality of the fits are presented by their x^/d.o.f. 
values. 



Coulomb correction factor 1/03(23) can be evaluated at every Qi, bin separately for the n n K 
and the n^n^n^ data samples. These correction factors are shown in Fig. ^ as a function of Q3 
in the range from 0.25 to 2.0 GeV. The resulting Coulomb correction applied to C^iQ^), due to the 
corrections to RsiQi) and Ri^iiQs), is shown in Fig. |^ where it is seen to rise as Qt, decreases, 
reaching the value of about 1 1% at ^3 = 0.25 GeV. 



5.3 The corrected 7i=^7i=^7i=^ BEC distributions 

The Coulomb corrected distributions i?3 (23), ^1,2(23) and Cti{Qt,) of the three-pion BEC, are shown 
in Fig. 0. A clear genuine three-pion Bose-Einstein enhancement is present in the low Qt, region, 
from about Qt, = 0.7 GeV reaching a value of C3(23) = 2.0 at Qj, = 0.2 GeV. The continuous line 
in the figure represents the fit result of C3(23) in the range 0.25 < Q3 < 2.0 GeV, as parametrised 
in Eq. |TT]. The quality of this fit, given by a %^/d.o.f. = 190/171, represents an improvement over 
the fit result obtained for the Coulomb uncorrected Ct,{Q3) distribution. The values obtained from 
the fit are = 0.580 + 0.004 fm for the emitter radius and a strength of X3 = 0.504 ±0.010, with 
correlation factor of p^.^ = +0.883. These are listed in Table [I] together with the fit results for K and 
£. A comparison between the results presented in the table shows, as expected, that the value of 
increases when the Coulomb correction is applied. We found that X3 increased by about 9% whereas 
r3 decreased by about 6%. 



5.4 Systematic errors 

To estimate the systematic errors we have considered the effects on the and results arising from 
the choice of the data selection criteria and from the procedure adopted for the Coulomb correction. 
We also investigated the effect on the results from our choice of the fitting range and the MC reference 
sample. These are summarised in Table ^ 

In order to estimate the systematic effects related to track and event selection the analysis was 
repeated restricting the track selection criteria described in Section 4.2. To evaluate the possible 
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Figure 3: (a) The Coulomb correction factors as a function of Q3 for the tt^tc^tc^ and for the tt^tt^ti^ 
distributions. The width of the bands corresponds to the statistical uncertainty, (b) The mean net 
Coulomb correction value as a function of Q2 for the genuine three-pion BEC distribution obtained 
by dividing the corrected C3{Q-i) distribution by the uncorrected one. 
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Figure 4: The Coulomb corrected tt^tt^tt^ BEC correlation distributions, (a) The non-subtracted 
^3(23) and 1,2 (23) distributions, and (b) the genuine Ct,{Q3) distribution. The solid line in (b) 
represents the fit results of Eq. [11] over the range 0.25 < Q3 < 2.0 GeV. 
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Fit variation 


^3 


Ah 


rs [fm] 


Ar3 [fm] 


a. The reference fit 


0.504±0.010 




0.580 ±0.004 




b. varying data selection cuts 

c. MC Coulomb correction 

d. fit range 0.25 < < 1-5 GeV 

e. fit range 0.30 < < 2.0 GeV 

f . addition of a long range term 

g. varying the MC reference sample 


0.519±0.017 

A /100_LAA11 

0.4oo ±0.01 1 

0.507 ±0.011 
0.471 ±0.011 
0.506 ±0.010 
0.496 ±0.026 


±0.015 

A A 1 
—0.016 

±0.003 

-0.033 

±0.002 

-0.008 


0.561 ±0.007 

A CO/; _L A AA/t 

0.586 ±0.004 
0.582 ±0.005 
0.569 ±0.005 
0.582 ±0.004 
0.597 ±0.015 


-0.019 

1 A AA/C 

±0.006 
±0.002 
-0.011 
±0.002 
±0.017 


Total systematic error 




0.041 




0.029 



Table 2: Results of several fits of the C2{Qt,) parametrisation to the data given with their statistical 
errors. The differences, AX3 and Ari,, between the parameter values of the reference fit (a) and the 
others from (b) to (g) are added in quadrature to obtain an estimate of the combined systematic 
uncertainty associated with the fitted and r3 values. These are given in the last row. 



contribution from our choice to account for the Coulomb effect on the BEC we also applied the 
Coulomb correction to the MC reference sample rather than to the data. Furthermore we considered 
for the fit two alternative ranges. We also investigated the influence of adding a quadratic 
term to the long range correlations. Finally, the systematic uncertainty coming from our choice of the 
JETSET 7.4 MC reference sample has been estimated by repeating the analysis with the HERWIG 5.8 



[ P8| ] generated sample. This alternative MC uses a totally different model of fragmentation (cluster 
fragmentation) from that used by JETSET 7.4 (string formation and fragmentation). An estimate of 
the over-all systematic uncertainties was obtained by summing in quadrature the differences between 
each fit (b)-(g) and the reference fit (a). 

The largest contributions to the over-all systematic error come from the choice of the selection 
criteria and from the choice of the lower fit range limit. As seen in Table ^ the results for r3 and 
X3 change by less than 4% when the Coulomb correction is applied to the Monte Carlo generated 
sample. We observe that is rather sensitive to the choice of the lower limit used as compared to 
the change of ri,. In addition to the list given in Table ^ we also investigated other possible sources of 
systematic effects, such as the choice of the Q3 bin size used in the fit, and verified that they indeed 
have negligible contributions. Finally, as noted above, the influence of the uncertainties of the ratio 
(n^^^/n^^^) on the subtraction formula and the fit results is also negligible. 

Thus the final values of the BEC parameters are: = 0.580 ± 0.004 (stat.) ± 0.029 (syst.) fm 
and X3 = 0.504 ±0.010 (stat.) ±0.041 (syst.), where the uncertainties of the measured parameters are 
strongly dominated by the systematic errors. In Fig. ^the 68% and 95% confidence level correlation 
contours for the BEC parameters are shown. The shape of the contours is determined mostly from the 
systematic errors. Accounting for the three-pion purity of 0.713 ±0.053, the BEC strength amounts 
to Xf'' = 0.707 ± 0.014 (stat.) ± 0.078 (syst.) for a 100% pure 7C±;c±7C± system. The purity error of 
0.053, due to the uncertainty of the MC generation rates, is incorporated in the systematic error of the 

pure 1 

A3 value. 
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Figure 5: The 68% and 95% confidence level correlation contours for the genuine three-pion BEC 
parameters and after Coulomb correction. The contours are calculated from the statistical errors 
and the systematic uncertainties listed in Table ||. The best values are represented by the solid circle. 

6 Relation to other experimental results 

A relation between the two-pion and three-pion emitter radii is derived in Ref . [Q] based on the Fourier 
transform of the source distribution which is assumed to be of a Gaussian shape. This relation confines 
the emitter range r3 as determined from a fit to 7?3(23): 

r2/V3 < ?3 < r2/V2 , (23) 

where is the two-boson BEC emitter size. Since in our analysis the emitter radius is determined 
from the genuine BEC distribution C^iQs), the previous bounds reduce to the equality: 

rs = r2/V2 . (24) 

In a former OPAL analysis [p3|], of the two-pion BEC present in the hadronic decays, two options 
were adopted for the reference sample. The first utilised the correlations of the pairs of tt+tt^ in 
the data and the second used the two-pion correlations of Monte Carlo generated sample. Since 
in our analysis the reference samples were taken from the Monte Carlo generated events, we have 
checked relation ( P^ with the previously measured OPAL value, obtained by the second method, of 
r2 = 0.793 ±0.015 fm, where only the statistical error was given. From Eq. ^it follows that this r2 
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value corresponds to rf-^'^ = 0.561 ±0.011 fm, which is in good agreement with our value obtained 
from a fit to the data of = 0.580 ± 0.004 (stat.) ± 0.029 (syst.) fm. 

The BEC of the n'^n^n^ system has also been studied by the DELPHI collaboration [ [131 ] using 
the hadronic Z° decays measured at LEP. In that analysis, which neglected the Coulomb correction, 
the following results were obtained: = 0.657 ± 0.039 (stat.) ± 0.032 (syst.) fm and = 0.28 ± 
0.05 (stat.) ± 0.07 (syst.). The relatively large statistical errors reflect the fact that a much smaller 
data sample was used than in the present measurement. Our result for r^, before Coulomb correction, 
of 0.616 fm with a statistical error of ±0.005 fm, is smaller than the DELPHI result but is consistent 
with it within one standard deviation. The measured X3 value depends on the pion track purity of the 
hadron sample analysed, so that caution has to be exercised when comparing the BEC strength values 
of different experiments. Keeping this in mind, we note that our mean X3 value lies considerably 
above that reported by DELPHI but it is still consistent with it within two standard deviations when 
the systematic errors are included. 



7 Summary and conclusions 

The Bose-Einstein correlations of three identical charged pions, produced in hadronic decays, 
have been studied after correcting for the Coulomb interaction. A significant genuine three-pion 
Bose-Einstein correlation signal is observed near threshold in the C3{Qt,) distribution obtained after 
the subtraction of the two-pion correlation contribution. The radius r^, of the three-pion emitter and 
the BEC strength are measured to be: 

rs = 0.580 ± 0.004 (stat.) ± 0.029 (syst.) fm and I3 = 0.504 ± 0.010 (stat.) ± 0.041 (syst.) 

where the uncertainties are dominated by the systematic errors. 

The Coulomb repulsive interaction opposes the Bose-Einstein enhancement in the low Q3 region 
and therefore it is reasonable that in our analysis the Coulomb correction increased the X3 value. This 
increase amounts to about 9%. On the other hand, the Coulomb correction has a smaller effect on 
the r3 value which is lowered by about 6%. Accounting for the three-pion purity the BEC strength 
amounts to ?if = 0.707 ±0.014 (stat.) ±0.078 (syst.) for a 100% pure n^n^n"^ system. 

A relation between the two-pion and the three-pion emitter dimensions was discussed in reference 
[^. We tested this relation by using the present result and that obtained in the latest OPAL two- 
pion BEC analysis [ [T5| ] where approximately the same data sample was used. The proposed relation 
between r2 and r^, expressed in Eq. is in good agreement with our results. 
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